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Abstract 

Supersymmetric scenarios incorporating thermal leptogenesis as the origin of 
the observed matter-antimatter asymmetry generically predict abundances of the 
primordial elements which are in conflict with observations. In this paper we pro- 
pose a simple way to circumvent this tension and accommodate naturally ther- 
mal leptogenesis and primordial nucleosynthesis. We postulate the existence of a 
light hidden sector, coupled very weakly to the Minimal Supersymmetric Standard 
Model, which opens up new decay channels for the next-to-lightest supersymmetric 
particle, thus diluting its abundance during nucleosynthesis. We present a general 
model-independent analysis of this mechanism as well as two concrete realizations, 
and describe the relevant cosmological and astrophysical bounds and implications 
for this dark matter scenario. Possible experimental signatures at colliders and in 
cosmic-ray observations are also discussed. 



1 Introduction 

Extending the Minimal Supersymmetric Standard Model (MSSM) with three heavy 
right-handed neutrino superfields is one of the best motivated scenarios for physics be- 
yond the Standard Model. The decoupling of the heavy degrees of freedom induces at low 
energies, after the electroweak symmetry breaking, tiny neutrino masses suppressed by 
the large right-handed neutrino masses; this is the renown see-saw mechanism Fur- 
thermore, the out of equilibrium decay of the lightest right-handed (s)neutrinos in the 
early Universe could have generated the observed matter-antimatter asymmetry through 
the mechanism of leptogenesis [2]. Successful thermal leptogenesis requires, though, a 
rather large mass scale for the new particles, M > 10^ GeV [3], which could destabilize 
the electroweak scale. Supersymmetry guarantees that the large quadratic quantum cor- 
rections to the Higgs mass introduced by the right-handed sneutrinos exactly cancel with 
the ones introduced by the right-handed neutrinos, thus avoiding the severe hierarchy 
problem of the non-supersymmetric version of the leptogenesis mechanism. 

It is remarkable that this simple scenario can simultaneously address two of the most 
severe limitations of the Standard Model, namely the existence of non- vanishing neutrino 
masses and the origin of the observed matter-antimatter asymmetry in our Universe (as 
well as the above mentioned hierarchy problem). This appealing scenario is nevertheless 
not exempt of problems. Gravitinos are very efficiently produced in the very hot plasma 
necessary to generate the observed baryon asymmetry through the mechanism of thermal 
leptogenesis, the relic density being 



where is the gravitino mass, rrig is the gluino mass and Tr is the reheating temper- 
ature of the Universe. If the gravitino is heavier than the lightest neutralino, it decays 
during or after Big Bang Nucleosynthesis (BBN). Being the gravitinos so abundant at 
the time of formation of the primordial elements, c/. Eq. (11.11) . the large hadronic energy 
injected into the primeval plasma destroys the successful predictions of the Standard 
BBN scenario unless the gravitino mass is larger than ~ 10^ GeV E] • 

On the other hand, if the gravitino is lighter than any observable supersymmetric 
particle, it constitutes a natural candidate for the cold dark matter of the Universe 
provided it is stable at cosmological scales. Namely, following Eq. (11.11) . the dark matter 
rehc density inferred by WMAP for the ACDM model, ^cdm^^ — 0.11 can be 
reproduced for the range of reheating temperatures required by thermal leptogenesis, 
Tr > 10^ GeV and for typical gluino masses, m-g ~ 1 TeV, provided the gravitino mass 
is larger than ~ 10 GeV. 
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This attractive scenario is in general in conflict with the observed abundances of 
primordial elements. If R-parity is conserved, the Lightest Observable Supersymmetric 
Particle (LOSP) can only decay into Standard Model particles and the gravitino with a 
decay rate suppressed by the Planck scale, the lifetime being: 



Therefore, the LOSP is typically present during or after BBN, jeopardizing the successful 
predictions of the standard nucleosynthesis scenario. This is in fact the case for the most 
likely candidates for the LOSP: the lightest neutralino and the right-handed stau. More 
precisely, when the LOSP is the neutralino, the hadrons produced in the neutralino 
decays typically dissociate the primordial elements [5l|8], yielding abundances in conflict 
with observations. On the other hand, when the LOSP is a charged particle, X~, the 
formation of the bound state (^HeX~) catalyzes the production of ^Li |9] leading to an 
abundance of ^Li in stark conflict with observations [TU] (for a recent review about BBN 
constraints see Ref. [TT]). 

Different solutions have been proposed to this problem. For instance, in some specific 
supersymmetric models the LOSP can be a sneutrino [12] or a stop [13], whose late 
decays do not substantially affect the predictions of BBN. For neutralino or stau LOSP 
a possible solution consists in introducing a small amount of i?-parity violation, so that 
the LOSP decays into two Standard Model particles before the onset of BBN, thus 
avoiding the BBN constraints altogether [2]. Maintaining the requirement of R-parity 
conservation, other solutions are to assume a large left-right mixing of the stau mass 
eigenstates [15], a LOSP mass that is nearly degenerate with the gravitino mass [16], or 
to assume some amount of entropy production after LOSP decoupling, which dilutes the 
LOSP abundance [IT]. 

In this paper we would like to propose a scenario which yields a thermal history of 
the Universe consistent with supersymmetric dark matter and with baryogenesis through 
thermal leptogenesis, without altering the successful predictions of the Standard BBN 
scenariolll We will assume the existence of a hidden sector fermion, X, lighter than the 
LOSP. Thus, new decay channels are possible for the LOSP, for instance, when the LOSP 
is the lightest stau or the lightest neutralino. 



If these decays are fast enough, the density of LOSPs at the time of nucleosynthesis 
can be significantly reduced and thus the successful predictions of the standard BBN 

^Sce also Ref. |18) for a related proposal with light axinos, however without supersymmetric dark 
matter. 
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scenario will not be jeopardized. If the gravitino is the LSP, the hidden sector fermion 
will eventually decay into the gravitino and other particles, hidden or observable. The 
late decays into gravitinos and hidden sector particles, if kinematically possible, may 
disrupt the abundances of hidden sector primordial nuclei, but not the abundances of the 
observed primordial nuclei. On the other hand, the decays into gravitinos and observable 
particles occur at a rate much larger than the age of the Universe, as we will show, 
thus not affecting primordial nucleosynthesis. Lastly, if the gravitino is not the LSP it 
will decay into hidden sector particles. Again, these decays may disrupt the primordial 
abundances of hidden sector nuclei, but will not have any impact on the standard BBN 
predictions. This mechanism is sketched in Fig.l, for the case where the gravitino is 
the LSP (left panel) and for the case where the gravitino can decay into hidden sector 
particles (right panel). 

The couplings of the hidden sector fermion to the MSSM particles are subject to a 
series of constraints which will be discussed in Section 2. In Section 3 we present two 
concrete models where the mechanism sketched above can be implemented. In Section 
4 we will comment on possible signatures of this scenario at colliders or at cosmic-ray 
observations. Lastly, in Section 5, we will present our conclusions. 

2 Hidden sector couplings to the MSSM 

We will consider in this paper a scenario where the MSSM particle content is extended 
with a light superfield (chiral or vector), which is a singlet under the Standard Model 
gauge group. We will further assume that the fermionic component of this superfield 
couples to the LOSP and its Standard Model counterpart via a tiny Yukawa coupling, 
hence we will refer to this fermion as "hidden fermion" . In this section we will carefully 
discuss the implications for leptogenesis of the existence of such hidden fermion, as well 
as the constraints on this scenario from BBN and from structure formation. Let us first 
discuss the case of stau LOSP and later on the case of neutralino LOSP. 

2.1 Stau LOSP 

The interaction Lagrangian between the hidden fermion and the lightest stau, ri, is given 
by the renormalizable term 

- C = Xr,XTTi+h.C. . (2.1) 
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Then the stau can decay either ri — )■ '(/'3/2 t or ri — > X r with decay rates: 

r~ . - ^ "^^^ (i "'/^V (22) 
^..x.-^^(l-4)^ (2.3) 

Sir \ J 

where the tau mass has been neglected. If the couphng A^i is large enough, the stau will 
decay before it can form bound states with ^He, thus preventing the catalytic production 
of ^Li. More concretely, the requirement that the lightest stau decays before ~ 2 x 10'^ s 
[HI [To] into Xt, implies the following lower bound on the coupling X-t-ti. 

If this condition is satisfied, the lightest stau will by itself not play any role during BBN. 

A necessary requirement for the viability of this mechanism is that cosmological 
constraints, namely from thermal overproduction of the hidden fermion X and from 
structure formation, are satisfied. Different scenarios can arise depending on whether 
the LSP is the gravitino or the hidden fermion, and on whether the NLSP is stable on 
cosmological time-scales or not. Below we will discuss each case separately. 



2.1.1 Cosmologically stable gravitino and hidden fermion 

If the hidden fermion is the lightest particle in the hidden sector, its decay channels 
into a gravitino and other hidden sector particles are kinematically forbidden (note that 
in this case the hidden fermion still can decay into a gravitino and Standard Model 
particles, e.g. X — )■ ip3/2T~^T~, X — ^^3/27, with lifetimes which can be larger than the 
age of the Universe, as will be discussed in Section 4). If this is the case, two particle 
species contribute to the dark matter, namely gravitinos and hidden fermions, each of 
them having a thermal component and a non-thermal component: 

^dm — ^3/2 ~l~ ^3/2 ~^ • (2-5) 

Here, ^3^2 non-thermal contributions to the gravitino and hidden fermion 

relic density, respectively, which are given by 

^3/2 = ^BR(n^^3/2r)n^^ (2.6) 
^5 = — BR(ri ^ Xt) nf^ , (2.7) 
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Figure 1: Sketch of our proposed mechanism, for the cases where the LSP is the gravitino 
{left panel) or a hidden fermion {right panel). The Lightest Observable Supersymmetric 
Particle (LOSP), the gravitino and the hidden fermion can be produced thermally or 
non-thermally through the decays of heavier particles; slow decays are indicated with 
dashed arrows and fast decays with solid arrows, whereas the other decay products 
are indicated with wiggled lines. If thermal leptogenesis is the correct mechanism to 
explain the observed matter-antimatter asymmetry, the LOSP decay into the gravitino 
occurs during or after the time of primordial nucleosynthesis, altering the predictions of 
the Standard BBN scenario. If this is the only decay channel of the LOSP the impact 
is usually dramatic, yielding abundances in conflict with observations. However, if the 
LOSP coupling to the hidden fermion is much larger than the coupling to the gravitino, 
this decay can occur before the onset of the nucleosynthesis reactions, thus avoiding 
altogether any possible effect of the LOSP on nucleosynthesis. Eventually the hidden 
fermion will decay into the gravitino and other hidden sector particles {left panel) or vice 
versa {right panel). Nevertheless, these decays do not alter the abundances of primordial 
elements in our observable sector. 

where is the stau thermal abundance 

~ 2 X 10^3 ( """^ ) ' , (2.8) 
VlOOGeV/ ' ^ ' 

corresponding to a yield Y^-^^ ~ 7x 10^^"^ (m^^^/lOO GeV) [I9]. Furthermore, 
tribution to the total dark matter density from thermally produced gravitinos, Eq. (11.11) . 
If the hidden fermion couples to the observable sector through a renormalizable coupling, 
the thermal production proceeds dominantly via the decay of thermally produced staus 
at temperatures T ~ rrir^. The corresponding hidden fermion relic abundance Q^x 
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given by [20 



^ 1.09xl0^'2m,vr,^.,v 

3/2 



100 \ / VTLx \ ( rn 



- 8.6 X 10^11. .)|A,|^^-J [^J[^-^J . (2.0) 

Here, (7* ~ 100 denotes the effective number of degrees of freedom at temperature T ~ 
rrir^ , and 6 parameterizes tlie potential enliancement of the hidden fermion abundance 
due to possible additional couplings of the hidden fermion to other MSSM particles, 
cf. Eq. fl3.9p below. Here we set for simplicity 6 = 0. 

In order to sufficiently reduce the number density of staus at the time of BBN it is 
necessary that BR(ri — )■ rX) ^ 1, and therefore, ^3^2 — 0- Requiring that the total dark 
matter density does not exceed the measured value by WMAP implies then 

+ + —^t ^ 0-11^^' • (2-10) 

As long as the thermal production of X and staus is small (which also implies a small 
coupling X^-^ ) this bound reduces to the standard overproduction constraint on thermally 
produced gravitinos, rig^g^^ ^ 0.11, which by its own allows high reheating temperatures 
Tr ^ 10^ GeV, as required by leptogenesis, provided the gravitino mass is large enough, 
see Eq. (11. ip . On the other hand, in the regime where the production of X is sizeable, 
Eqs. (12. 9p and (I2.10p combine to a bound on the coupling 

|A,l<10--(^)"7l-^)". (2,11) 

\mxj V rnij 

which is independent of the gravitino mass. 

The results are illustrated by the red lines in Fig. [2] for the case of gravitino LSPs 
and in Fig. |3]for the case of hidden fermion LSPs. The lines show, for different reheating 
temperatures and different masses of the hidden fermion, the value of the coupling A^^ 
as a function of the gravitino mass from the requirement that the total dark matter 
density is equal to the value inferred by the WMAP collaboration. In these plots, it 
was assumed for definiteness m^^^ = 250 GeV. On the other hand, the shaded regions 
correspond to choices of parameters where the stau decays with a lifetime longer than 
2 X 10^ s, thus leading to ^Li overproduction. As apparent from the plots, in both cases 
of hidden fermion LSP and gravitino LSP, there is a fairly wide region of the parameter 
space where the reheating temperature can be large enough to allow thermal leptogenesis 
while preserving the successful predictions of the standard BBN scenario. 
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Figure 2: Summary of constraints on stau LOSP scenario with gravitino LSP, as de- 
rived from Eq. (I2.ip . as function of the stau-tau-hidden fermion Yukawa-couphng A^^ 
and gravitino mass 7713/2. The masses of X and ri are fixed as indicated. The red and 
blue lines show the values of ^3/2 and A^^ that yield the correct total rehc abundance 
for different reheating temperatures Tr = 10^-10^ GeV, assuming that X is stable or un- 
stable, respectively (the gluino mass has been set to rrig = 800 GeV). Furthermore, the 
dashed part of the blue lines is excluded by constraints on mixed warm/cold dark matter 
as discussed in the text, and the gray region is excluded by ^Li overproduction during 
BBN. It is clear that for vanishing coupling to the hidden sector, reheating temperatures 
around 10^ GeV, as required by thermal leptogenesis, are in conflict with BBN whereas 
an allowed window opens up for non-zero A^^. 



2.1.2 Unstable hidden fermion 

If the gravitino is the LSP, and if kinematically allowed, the hidden fermion X decays 
into gravitinos and hidden sector particles well before matter-radiation equality. Note 
that the particles produced in the decay interact very weakly with the particles in the 
observable sector, therefore the late decays of the hidden fermion do not modify the 
abundances of primordial elements in the observable sector. In this scenario, the dark 
matter consists of thermally produced gravitinos, with a relic density given by Eq. (11.11) . 
and non-thermally produced gravitinos, coming from the late decay of hidden fermions 
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X and staus ri. The dark matter abundance is then given by 

^^dm = n'i), + nl)^ + ^ (n'^ + fil) , (2.12) 

=f2-VVDM 

where we assumed for simphcity that the hidden-sector particles produced in the decay of 
X are massless or very hght and hence contribute neghgibly to the rehc abundance|§ The 
component coming from the late decay of X, as well as the small fraction of gravitinos 
produced directly in ri decays, will typically act as warm dark matter (WDM), with free- 
streaming lengths Aps ^ 5Mpc (c/. Fig. 4 in Ref. [2T])|f| Here, the free-streaming length 
of a particle is defined as the distance the particle has traveled between its production 
and the onset of structure formation. It is given by 

Zp 

Afs= I d^^y (2-13) 

3000 

where Zp denotes the red-shift at which the particle is produced and is a function of the 
lifetime of the parent particle, z ~ 3000 is the redshift at matter-radiation equality, v{z) 
is the particle's velocity and H{z) denotes the Hubble parameter as function of redshift, 
see e.g. Ref. [22] . 

Observations of the power spectrum of high-redshift Hydrogen clouds via the Lyman- 
a forest [23] imply the upper bound Aps ^ 0.5 Mpc [21], when all dark matter components 
have comparable free-streaming lengths. This bound relaxes if a large fraction of the dark 
matter is cold and just a small fraction of it is warm. Bounds on the fraction / of the dark 
matter density that is allowed to be warm with a free-streaming length above 0.5 Mpc, 
were discussed in Ref. [25] in the context of sterile neutrinos. There, using Lyman-a 
data [26] and WMAP5 results, 2cr-bounds around / < 0.05 were found for a warm 
component with free-streaming lengths around 0(10 Mpc), corresponding to 0{1 km/ s) 
thermal velocities. 

For definiteness we will take / = 0.05 throughout this paper, which implies, following 
Eq. ( ]2.12p . the requirement of a small thermal abundance of ri and X. Allowing a fraction 
/ of dark matter to be warm, and provided that ^3^2 — O5 gives then 

^fe + f^S) </0.11/^-^ (2.14) 



^Otherwise they would contribute to the warm dark matter component and they could cause dan- 
gerous late decays into Standard Model particles. 

■^Note that this could also be relevant for J in the above case where gravitinos and hidden fermions 
are stable, since for small hidden gaugino masses the free-streaming length becomes large. However, in 
the example shown in Fig.[2l the red lines feature always free-streaming lengths below Aps ^ 0.5 Mpc in 
the region where A^i ^ 10"^"*. Furthermore, since the stau yield is small, the impact of this component 
on mixed warm/cold dark matter bounds is negligible even if the free-streaming lengths are large. 
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Figure 3: Summary of constraints on stau LOSP scenario with hidden fermion LSP and 
gravitino NLSP, similar to Fig. El Here, the mass of the hidden fermion X is assumed 
to he in the range mx = 0.1-10 GeV, whereas the reheating temperature and the ri 
mass are fixed as indicated. The red and blue lines show the constraints for the case of 
a stable or unstable gravitino, respectively. Dashed parts of the lines are again excluded 
by constraints on mixed warm/cold dark matter. 

In addition to this constraint, the viability of the present scenario requires that both, the 
thermal gravitino abundance and the thermal hidden fermion abundance, do not exceed 
the total dark matter density, namely Vt^^j^ + (m3/2/mx)f^x ~ O-^l ^/ ^q- (12. 9p . 

The impact of all these constraints on the parameter space of the scenario is shown 
in Fig. |2] as blue lines. The dashed part of the lines is excluded from the constraints on 
mixed warm/cold dark matter, Eq. f l2.14p . Compared to the scenario where the hidden 
fermion is stable, now larger values of the coupling A^^ are allowed, since the hidden 
fermion does not directly contribute to the total dark matter abundance any more. 

2.1.3 Unstable gravitino 

In some scenarios the hidden fermion could be the LSP. If this is the case, the gravitino 
can decay into it, if kinematically allowed, yielding a scenario which is qualitatively 
different to the one studied in the previous subsection. The total dark matter abundance 
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is given in this case by: 

^dm = + "I ( ^3/2 + ^3/2 ) ' (2-15) 

m3/2 V ^ 

^ V ' 

=f^WDM 

where we assumed that hidden-sector by-products of the gravitino decaying into X are 
masslessEI In this case, the thermal production of gravitinos itself produces ultimately 
WDM, yielding the constraint 

^t/2h' < ^ / 0.11 , (2.16) 

where we again assumed that ^3^2 — 0. It is apparent from this equation that in scenarios 
with unstable gravitinos and high reheating temperatures, small masses mx ^ / ^123/2 
are favored, to avoid strengthening the bounds on thermal gravitino production. Fur- 
thermore, the thermal or non-thermal production of X becomes essential, since it must 
yield the dominant part of the dark matter abundance according to + ^ h^'^ 0.11, 
which implies a coupling to the hidden fermion like Xr^ ~ 10~^^(m?^j/mj(:)^^^, as long as 
fij is negligible. 

This situation is illustrated by the blue lines in Fig. [3l for fixed reheating temperature 
Tr = 10^ GeV and for different masses of the hidden fermion mx- As apparent from 
this plot, for being in agreement with the bounds on mixed warm/cold dark matter 
the gravitino mass has to be considerably higher than mx in each of the shown cases, 
e.g. for mx = 1 GeV the gravitino mass should exceed 20 GeV. Typical couplings where 
the mechanism works lie in the range A^^^ ~ 10~^^-10~^*^. 



2.2 Neutralino LOSP 

We will now briefly discuss the case of neutralino LOSP. For definiteness, we consider an 
effective interaction Lagrangian between the hidden fermion X and the lightest neutralino 
Xi given by 

-C = gHXx',h' + gzX 7^x? Z^^ + g,X j,Xi A^^+h.c. . (2.17) 

The couplings gi include gauge couplings, weak mixing angles etc. The neutralino can 
decay either Xi ^ (^°,7, /i°)V^3/2 or xi ^ (^°,7, h'^)X. 

The electromagnetic and hadronic energy eyis = ("^^0 — ^i>3/2/x + "^l//i)/(2"^xi') 
released during these decays can induce photo- and hadrodissociation of ^He during 

"^Note that here again the component has a large free-streaming length Ai?5 > 1 Mpc in some 
cases, depending on the coupling A^i and the mass 111^^ . However, due to the small stau (and later the 
neutralino) yield, the impact of this component is always negligible in our plots. 
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BBN, and can lead to inter-conversions of protons and neutrons [Sj |8] . Without hidden 
fermion X (corresponding to the hmits Qi 0), this leads to stringent constraints on the 
gravitino mass and consequently also on the reheating temperature [27] . However, if the 
couplings Qi are large enough, the neutralino decays into hidden fermions already before 
BBN, allowing a high reheating temperature Tr ~ 10^ GeV as required by leptogenesis. 

As an example we will consider the decay Xi ~^ Z^X, where the decay width is given 

by 



9% 
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(2.18) 



assuming that only gz 0. The boson then decays into hadrons with branching ratio 
Bh ~ 0.7. 

For a reference neutralino yield of e^jj^y^o ~ 100 GeV X 10-^2 [28|, the most stringent 
constraints come from overproduction of ^He due to interconversion processes, as well 
as D production by hadrodissociation, leading to an upper bound on the neutralino life- 
time in the range r^o < Tmax ~ 1 — 100 s O [8]. In the parameter region consistent with 
thermal leptogenesis, the decay into gravitinos is negligible, and BR(x'j' — ?■ Z^X) 1. 
This implies a lower bound on the coupling, which reads in the limit m^o ^ mz,mx'- 



Xi 



200 GeV 



-3/2 



VlOOs 



-1/2 



(2.19) 



constraints for K,o 

Al 

decay modes Xi ~^SZ^,'~f)'ip3/2 into account using the rates given in [2 



The constraints from overclosure and from free-streaming are very similar to the 
case of stau LOSP discussed before, and we do not repeat them here. The corresponding 

10~^^ are summarized in Fig. HI In this figure, we also took the 

assuming a 

bino-like neutralinojfl Note that for both cases, gravitino and hidden fermion LSPs, one 
can find parameters compatible with Tr ~ 10^ GeV, and typical mixing parameters lie 
in the range gz ~ 10~^'^-10~^'^. 

The bounds shown in Fig. H] are relying on the relatively small adopted neutralino 
yield, and they can change qualitatively when the yield is much larger. For example, in 
the case where the hidden fermion is the NLSP and unstable (blue lines in the left panel of 
Fig. H]), the warm dark matter component produced by the decay chain Xi ~^ ^ ^ ^3/2 



^The BBN bounds were obtained by interpolating between the bounds for the different hadronic 
branching ratios given in [5]. 



12 



10-" 



= 50 GeV ni„e„Mi.„ = 200 Ge V j^„„ = 1 0" 




Tj = 10' GeV 



10" 10' 
[GeV] 



r„ = 10' GeV "i„™,„ii„„ = 200 GeV 



10"" - 



10"' 



unstable eravilino 



OA GeV 




10" 10' 

m„2 [GeV] 



Figure 4: Summary of constraints on neutralino LOSP scenario with gravitino LSP [left 
panel, as in Fig. |2]), and hidden fermion LSP [right panel, as in Fig. [3]), as function of 
the Xi^^^ couphng gz and gravitino mass 7713/2, and for Y^o = 10~^^. The constraints 
from BBN obtained in Ref. [S] [shaded region) exclude neutrahno hfe-times r^o > 100 s. 
However, Ref. [5] presents somewhat stronger bounds, of the order 1 — 10 s [grey dashed 
lines). 



imphes an upper hmit on the neutrahno yield. 



-12 20 GeV / 



m3/2 0.05 



(2.20) 



which follows from the bounds on mixed warm/cold dark matter. On the other hand, in 
the cases where the gravitino is the NLSP, one can find viable scenarios even for much 
higher values of the neutralino yield, which however requires that the hidden fermion is 
very light. 



3 Examples 

After having discussed the proposed mechanism as well as astrophysical and cosmological 
bounds in the last section in general, we will now present two concrete scenarios with 
hidden fermions from a vector and from a chiral supermultiplet, respectively. 
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3.1 Hidden gauginos of an unbroken [7(1) x 



Let us consider the case where the hidden fermion arises from the gaugino component of 
a vector superfield of an unbroken hidden U{l)x symmetry, which mixes with the U{1)y 
of hypercharge via a small kinetic mixing x <^ 1 (for details see also Ref. [2T]). This 
scenario is an example for the case where the superpartner of the hidden fermion X, 
here the U{l)x vector boson, remains exactly massless. For simplicity, we will assume 
that all matter charged under U{l)x is vector-like and heavy enough to be cosmologically 
irrelevant. 

When SUSY is exact, the canonical normalization of the kinetic terms of the vector 
superfields produces an unobservable shift of the hypercharge gauge coupling, while the 
hidden U{l)x gauge boson and gaugino completely decouple from the observable sec- 
tor [29]. However, in presence of SUSY breaking effects, the decoupling of the gaugino 
is not complete any more [21]. In the component formalism, the relevant part of the 
Lagrangian, including the supersymmetry breaking soft masses, reads 




= -^[X,. B^,) >C - , [\x \B)IC^"d,\ (3.1) 

where /C and Ai denote, respectively, the kinetic and mass mixing matrices 

IC=(' A and M=(^- , (3.2) 

and Xx/B and Dx/b are the gauginos and D-terms corresponding to the gauge fields Xfj,u 
and B^^, respectively. Note that the generation of mass mixing in general depends on 
details of the underlying theory. 

In the basis where the kinetic terms are canonical, which can be achieved by the 
redefinition (which also holds for the corresponding D-terms and the vector bosons) 





(3.3) 

the hidden gaugino and the four MSSM neutralinos mix, and the corresponding extended 
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(5 X 5) neutralino mass matrix reads, to lowest order in x 

(Mx 
6M 




5M 
Mb 





V 



-MzCfjSw 
MzSpSw 






Mw 
Mzcpcw 
-MzSpCw 





-MzCpSw 
MzCpCw 




\ 

MzSpSw 
-MzSpCw 

/ 



(3.4) 



where 5M ~ 5M-xMx, Mx ^ Mx and Mb ~ Mb- Here, /i denotes the MSSM /i-term, 
the mass of the gauge boson, sy/ the sine of the Weinberg angle and is related 
to the ratio of the two Higgs VEVs. 

Upon diagonalization of the previous mass matrix we find a mass eigenstate X = 
Xx + 6Ab with mass mx — Mx, namely it is mostly hidden fermion with a small bino 
admixture given by 

5M 



e 



Mb - Mx 



(3.5) 



from which follows that typically O ~ 0{x)- Thus, even though the canonically nor- 
malized fields X^^ and Dx completely decouple from the observable sector, the mass 
eigenstate X, couples to the observable sector {e.g. the right-handed stau tb) through 
the tiny bino component 



C D V2g%^i\BPRr)rR ^ V2g%^eiXPRr)rR . 



(3.6) 



The presence of an additional U{l)x gauge group opens the stau decay channel 
Ti — tX (we assume that ti ^ tr), with lifetime 



ge^y^™^ 1 1 



m 



X 



2 

n 



(3.7) 



In contrast to the scenario discussed in Sec. 2, where only the minimal effects of the 
coupling (12. ip where taken into account, many additional channels for thermal production 
are now open, since the U{l)x inherits all couplings of the hypercharge bino. The thermal 
production of hidden gauginos from 2 — 2 scattering processes was calculated in Ref. [21] 
and yielded an upper bound on the mixing angle G that is given by (provided that X is 
stable) 



e < 3 X 10 



-12 



ma 



mx 



(3.8) 



where mq denotes the squark masses. In this calculation the dominant scattering pro- 
cesses involving one QCD and one hypercharge vertex were taken into account. Addi- 
tional contributions from 1—7-2 processes were discussed in Ref. [20], yielding a relic 
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density according to Eq. (12. 9p . when identifying \~ 



= 9' 



'2 Yi and 




) 



2 



Y? m? \ 1 — 171% I m~ 



(3.9) 



where / runs over all sfermions of the MSSM. For a typical SUSY mass spectrum, one 
obtains 5 ~ (9(3-10). For stable X, this yields an upper bound of 



V fnx 

Note that both contributions from 1 — 2 and 2 — 2 processes are roughly of the same 
order of magnitude and not much stronger than the production coming only from a term 
like Eq. (O) alone. 

Bounds on the mixing parameter G are illustrated in Fig. [5l From there it is clear 
that a kinetically mixed hidden U (l)x gauge group with mixing parameters 10"^'^ ^ © ^ 
10~^° satisfies all constraints. Apart from the somewhat stronger thermal production, the 
situation is similar to what is shown in Fig. |2] and [31 

The required mixing parameter ~ x lies in the broad range of values that can be 
accommodated in string motivated \J{\)x extensions: Without additional symmetries, 
the kinetic mixing x is generically generated on one-loop level by integrating out chi- 
ral superfields charged under both, visible and hidden sector. In this case it acquires 
values typically around x ~ 10~'^-10"2, corresponding to one-loop suppression [291 130] . 
However, in compactifications of heterotic and type II strings, much smaller mixings are 
possible [311 |32l [331 [311 135] . For example, a lower bound around x ^ 10~^^ was argued to 
hold in cases of gauge mediated supersymmetry breaking in heterotic string models [3T] . 
whereas in type-II models with warped extra dimensions the kinetic mixing parameter 
could be parametrically even smaller [32]. Hence, scenarios with hidden unbroken \J{\)x 
gauge groups provide simple and natural scenarios where the tension between thermal 
leptogenesis and gravitino dark matter is solved [2T] . 

3.2 A model for small couplings from non-renormalizable op- 



Let us now consider the case where the hidden fermion X is part of a chiral supermultiplet, 
singlet under the Standard Model gauge group. Then, the gauge symmetry allows the 
superpotential term 




(3.10) 



erators 




(3.11) 
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Figure 5: Constraints on the parameter space of the hidden U{l)x gaugino scenario with 
a Ti LOSP, for the case of gravitino LSPs {left panel, of. Fig. [2]) and for hidden gaugino 
LSPs {right panel, of. Fig. [3]), as function of the mixing parameter 6 and the gravitino 
mass m3/2, and for ti mass fixed as indicated. In the left panel the reheating temperature 
varies and whereas in the right panel the hidden gaugino mass mx varies. In both cases 
the NLSP is unstable and decays into the LSP, since the hidden vector boson remains 
massless. 

This term leads, after the electroweak symmetry breaking, to a Lagrangian term of the 
form Eq. (12. ip . with a Yukawa coupling which reads 

An~^^~10 «33Cos/3^^^^J . (3.12) 

It is remarkable that this value for the coupling lies in the allowed region of Figs. [2] and |3] 
if the non-renormalizable operator is suppressed by masses close to the Grand Unification 
scale. 

The superpotential Eq. (l3.111) can naturally arise from the decoupling of heavy par- 
ticles with a mass of 0{M). Let us consider an extension of the MSSM with three 
additional chiral superfields: X (1,1,0), H'^ (1,2,1) and H'^ (1,2,-1), where in paren- 
thesis we indicate the quantum numbers under SU{3)c x SU{2)l x U{1)y- Note that H'^ 
and H'^ have identical gauge quantum numbers as the MSSM Higgs doublets Hu and H^, 
respectively. To avoid unwanted terms, we will further impose the following Peccei-Quinn 
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transformation on the superfields: 



{Q,U^^,D^^,L,E%,m^) ^ e^-{Q,U^^,D%,L,E%,N^^), (3.13) 
{H^,H,) ^ e"2-(i/„,i7,), (3.14) 

K ^ e^'^'K, (3.15) 
H', ^ e-2^"i/^, (3.16) 
X ^ X. (3.17) 

This symmetry forbids the bihnear term H^Hd in the superpotential as well as a Majo- 
rana mass term for the right-handed neutrinos. The bilinear fi term could be generated 
via the Giudice-Masiero mechanism [36]. On the other hand, in order to generate the 
right-handed neutrino masses, we will further introduce a Standard Model singlet, $, 
which transforms under the Peccei-Quinn symmetry as $ — > e~^*°$. Then, the term 
^N^N^ is allowed in the superpotential and leads to right-handed Majorana masses if 
the scalar component of $ acquires a vacuum expectation value, ($) ~ M/j. Note that 
even with the presence of the new field it is not possible to generate a bilinear term 

Then, the renormalizable superpotential of the model reads 

W = L{Xe ■ Hd)E% + Q{Xd ■ Hd)D'j, + \uQHuU% + KLHuN'j, + Xm^N^N''^ 

+ {^JiD ■ Hd)K + (Ax ■ Hd)KX + MxXX + KxX^ + CxX + . . . , (3.18) 

where the ellipsis indicates additional terms in the superpotential which lead to the 
breaking of supersymmetry and the breaking of the Peccei-Quinn symmetry. Here, we 
used a compact notation where iJ^i = (if^, H'^, Xe-Hd = XeHd+X'^H'^, Xx-Hd = XxHci+ 
X'^H'^ etc. The linear term in X can be eliminated by a shift, X — X — Cx/(2Mx). Upon 
redefining the parameters, this amounts to setting Cx = 0. In general ~ C'(^gut), 
although we will choose to work in the basis in the {Hd-, -ff^)-space such that /i^i = (0, M), 
being M ~ (9(Mgut)- Written in components, we thus have 

W = W^^^^ + X,LHuN% + Xm^N%N% 

+ X',LH',E% + X'dQH'dD'n + MH'.H'^ + (A.i/, + X'MKX 

+ MxXX + KxX'^ + ... . (3.19) 



®In the presence of the field 4> new, non- renormalizable, terms appear in the superpotential suppressed 
by a large mass scale, M*, such as QQQL^'^/M^. After the breaking of the Peccei-Quinn symmetry 
this term leads to a superpotential term {Mii/M^)'^QQQL/M^ which induces proton decay. However, 
the small factor {Ma/M-tY ^.nd the plausibly small coefficients of the dimension-7 operator for the first 
generation can yield a proton lifetime in agreement with the stringent experimental bounds. 
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The breaking of the Peccei-Quinn symmetry at intermediate scales leads to Majorana 
masses for the right-handed neutrinos, Mm- Then, the decoupling of the heavy fields H'^, 
H'^, finally leads to the following effective superpotential 

WeE =WS^^ - {KM]^^\l){LH,){LHa) 

- ^LH^E^X - ^QH^DIX + MxXX + kxX^ , (3.20) 

which contains the term Eq. (13.111) . Note that this model is free from gauge anomalies 
and preserves the successful MSSM gauge coupling unification. 

In this model, there are additional contributions to the thermal production of the 
hidden fermion X from 2 — )■ 2 scattering processes, such as r /i*^ — X r and t t ^ X hP. 
These processes stem from the higher-dimensional operators (HDO) in Eq. ( 13.20p . and 
the amplitude for these processes increases with energy squared, IM-Wdq ~ s/M"^, where 
s denotes the square of the center-of-mass energy. Therefore these production channels 
can be very efficient at the high temperatures required by thermal leptogenesis. We 
will now shortly discuss the impact of these production channels on the overproduction 
constraints. 

The sum of the matrix elements for all possible scatterings producing an X or an X, 
respectively, are 

channels channels i,j 

where we also summed over the initial and final-state spins. The effective X-stau-tau 
coupling can be identified with 

^?.-#E(^^)i-^. (3-22) 

i 

where = v ■ cos /3 ~ cos /3 175 GeV. Assuming for simplicity that all other Yukawa cou- 
plings are zero (which will give the minimal contribution to the UV thermal production), 
the hidden fermion yield reads [20] 



where Mp ~ 2.4 x 10^^ GeV denotes the reduced Planck mass. In this expression the 
prefactor 3 takes into account that the hidden fermions are produced in two ways: the 
scattering 12 — j- 3X, which produces just one hidden fermion, and the scattering 12 — t- 
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3X, which produces two hidden fermions, due to the fast decay X XX (unless kx ^ 
10^^^). Assuming that X is stable on cosmological time-scales, this leads to an abundance 

.™.^.0W(^)(^)(15^)^(^)-. (3.., 

For comparison with the freeze- in abundance (12.9p . it is instructive to express this in 
terms of the effective coupling parameter A^r^. Then, with = 915/4, 

.-.-o,nco.-/,(^) (^) [^^y . (3.5) 

From requiring < 0.11 one thus obtains a bound 

|A,|<2.7xlO-"cos,(i^)'(15^)\ (3.6) 

which depends on the reheating temperature Tr and the mass of the hidden fermion mx- 

Bounds on the above model, for the case where the hidden fermion mass is very small, 
are summarized in Fig. [6] (note that gravitino LSPs together with reheating temperatures 
Tr ~ 10^ GeV are excluded in the present model due to overproduction of the hidden 
fermion NLSPs). We show results for unstable gravitinos only, since even a very small kx 
makes the gravitino unstable at cosmological timescales, due to the fast decay ?/'3/2 — > 
XXX. When comparing Fig. [6] with Fig. |3l it is apparent that the contribution to the 
relic abundance of hidden fermions coming from the HDO, cf. Eq. (13. 23 p . reduces the 
allowed parameter space for the coupling A^^^ considerably. These ultraviolet contributions 
depend linearly on the reheating temperature and dominate the infrared contributions 
coming from the renormalizable operator. They can potentially reintroduce the gravitino 
problem, cf. Eq. (II. ip . if the mass of the hidden fermion is too large. However, for small 
hidden fermion masses nix ^ 1 GeV all constraints from overproduction, BBN and 
structure formation can be simultaneously satisfied for couplings in the range ~ 
10^^'^-10~^'^. Most interestingly, these couplings are preferred if the non-renormalizable 
operator is suppressed by masses close to the GUT scale, see Eq. (I3.12p . 

4 Experimental signatures 

In this section we will briefly discuss possible experimental signatures of the scenario 
proposed above at colliders and in cosmic-ray observations. 

In this scenario the coupling constant of the LOSP to the hidden fermion X is 
typically smaller than 10~^^, therefore the decay length is much larger than the size of 
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Figure 6: Similar to Fig.[3l but for the model described in Section. [3l2| taking into account 
also the minimal ultraviolet contributions to the X production, and assuming tan /3 = 2. 
We show constraints for the case where the gravitino is the NLSP and decays into X 
(in contrast to the previous figures the gravitino decays here into three hidden fermions, 
^3/2 XXX). The reheating temperature and the Ti mass are fixed to 10^ GeV and 
250 GeV, respectively, whereas the hidden fermion X mass varies in the range mx = 
0.01-10 GeV. Note that a typical value for A^i in our model is ~ 10~^^-10~^^ when 
the non-renormalizable operator is suppressed by masses close to the GUT scale, cf. 
Eq. (l3:T2|) . 



typical collider detectors. As a consequence, if the LOSP is the lightest neutralino, the 
experimental signatures at colliders are identical to the case of the MSSM with R-parity 
conserved, since the lightest neutralino just escapes the detector. On the other hand, 
if the LOSP is the lightest stau, it propagates through the detector leaving a heavily 
ionizing charged track [37] . Furthermore, if the stau velocity is small enough, it could get 
trapped in the detector and decay eventually, producing a tau moving in a non-radial 
direction. This signature, albeit very spectacular, is not specific of this scenario but also 
arises in scenarios with stau LOSP and gravitino or axino LSP [38] . 

A hint towards our proposed scenario arises from the measurement of the coupling 
stau-tau-hidden fermion. More concretely, at colliders it will be possible a determination 
of the stau mass and the stau lifetime, from which the coupling could be determined 
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through Eq. fl23|) : 



If the couphng inferred from colhder experiments hes in the range preferred by cosmology, 
Eq. (12. lip , and if the hfetime Tfj is in agreement with BBN bounds, this scenario would 
gain strength. 

Furthermore, if the small coupling of the stau to the tau and the hidden fermion 
is due to the non-renormalizable operator Eq. (l3.1ip . then a very interesting signature 
could be observed at colliders: after being stopped in the detector, the stau could decay 
producing also a Higgs boson, with branching ratio: 



, ^, tan'' (3 sin a m~ o/ 2^ r a 

BR(ri ^ ThX) ^ ——-^-±f{ml/ml) , (4.2) 



where /(r) = 1 + 9r — 9r^ — + 6r(l + r) ln(r), v = 175 GeV is the Higgs vacuum 
expectation value and a is the angle that diagonalizes the CP-even Higgs squared-mass 
matrix. For typical values of the light Higgs mass and the stau mass, we obtain 

BR(?.^.«). 0.006 (^)'(^^)\.n^., (4.3) 

where we have taken mh/mr^ = 115/250 in the argument of the function /. Therefore, 
if 1000 staus could be stopped, a few Higgs events could be observed in the detector. 

For comparison, in scenarios with gravitino LSP and stau NLSP, the Higgs decay 
is strongly suppressed by the tau Yukawa coupling. The branching ratio, taking Higgs- 
strahlung into account, reads 

, tan^ B sin^ a mi ,9,0, , , 

BR(ri ^ r/1^3/2) ^ -^g{ml/ml) , (4.4) 

where g{r) = (1 + 8r + 6r^) ln(r~^) — (43 + 80r — 108r^ — 16r^ + r^) /12, and we assumed 
m3/2 ^ m^j:-^^. Again, inserting typical values of parameters, 

/tan /3\ ^ 

BR(?i ^ r/i^3/2) ~ 3.4 X 10"' f \ sin^ a , (4.5) 

which is suppressed by the tau mass, making the observation of Higgs events from stopped 
staus more difficult. Thus, the number of Higgs events in late stau decays constitutes a 
way to discriminate between a scenario with gravitino/axino LSP and a scenario where 
a hidden fermion couples to the stau and the tau via the non-renormalizable operator 
Eq. dHUD. 



Another possible experimental signature in this scenario is the observation of a large 
number of lepton flavour violating stau-LOSP decays. In the model described in section 
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13.21 it is in general not possible to simultaneously diagonalize the usual MSSM cou- 
pling XgLHdE'^ and the new couphng with the heavy Higgs doublet X'^LH'^E]^. As a 
consequence, the effective coupling in the Lagrangian {\'^ij\x{Hd) / MLiE'^^X in gen- 
eral contains sizeable off-diagonal entries, which induce at tree level the lepton ffavour 
violating decays Ti — ?■ and Ti — )■ eX with branching ratio: 

BR(ri ^ IjX) 
BR(?i ^ rX) 

Future colliders are capable of detecting the electrons and muons produced in the decay of 
long lived staus if (A;;)i3/(A;,)33 is larger than ~ 3x 10 (9 x IQ-^) provided 3 x 10^ (3x10^) 
staus can be collected [32] • Therefore, this scenario offers good prospects to detect flavour 
violation in stau decays at future colliders. 

It is important to note that this scenario can be easily compatible with the present 
experimental constraints on lepton flavour violation in the charged lepton sector. Namely, 
the process /i — )■ 67 induced at the one loop level via the Yukawa coupling Ag is strongly 
suppressed due to the large mass of H'^. Furthermore, if the scale of mediation of SUSY 
breaking is larger than the mass of the heavy Higgs doublets, which is the case for 
^^3/2 ^ 10 GeV(M/10^^ GeV), the couplings Ag and A^ induce off-diagonal entries in 
the squark and slepton mass matrices through renormalization group running, which in 
turn induce flavour violating processes in the quark and lepton sector through quantum 
effects. However, these effects appear then at the two loop level and are naturally sup- 
pressed. The most stringent bound coming from // — )■ 67 implies mild constraints of the 
order (Ag^Ag)2i < 0.1 (m^:^/250 GeV)^, whereas the 31 and 32 entries can be order one. 
Moreover, these effects can be further suppressed by appropriate choices of the flavour 
structure of the matrices Ag and A^. Therefore, in the scenario proposed in section 13.21 
it is expected that the stopped staus will decay not only into taus, but also into muons 
and electrons with a large rate, while being consistent with all present constraints on 
lepton flavour violation. In contrast, in scenarios with gravitino LSP and stau NLSP the 
rates of lepton flavour violating stau decays are predicted to be small [39]. Therefore, 
the number of lepton flavor violating stau decays offers a sensitive probe to discriminate 
gravitino LSP scenarios from the hidden fermion scenario discussed in section 13. 2[ 

In section [2.1.1l we discussed the possibility that the gravitino and the hidden fermion 
could be stable at cosmological timescales. Nevertheless, the heavier dark matter particle 
will eventually decay into the lightest and a tau-antitau pair with a decay rate which is 
doubly suppressed by the Planck mass and by the small coupling between the stau and 
the hidden fermion. Namely, when the gravitino is lighter than the hidden fermion, the 



i3 



33 



(4.6) 
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latter decays with a rate which reads 

r(X = ^^^gj;3M^g^,„^ . (4.7) 

where we have assumed 7713/2 <^ ttlx- On the other hand, when the hidden fermion is 
hghter than the gravitino, the decay rate approximately reads [H 



I \~ I ^ Tn7 

rW,3/.-r-r-A-).^^. (4,8) 

The electrons, positrons, gamma rays and neutrinos produced in the decay could be in 
principle detected in cosmic-ray observations as an excess over the expected astrophysical 
backgrounds. For dark matter particles with a mass ~ 100 GeV, the decay products are 
observable if the lifetime is shorter than ~ 10^^ s |10]. However, in our scenarios, the 
heavier component of dark matter has a lifetime 

.(X-..V-.3.)^I0».(J^)-^(^)-'(^)^(^)^ , ,«) 

^.-A-)^:03».(J^)''(^)-'(^)* , ,4.0) 

which have been normalized to typical values of the masses and couplings which yield a 
thermal history of the Universe incorporating thermal leptogenesis and successful BBN 
(c/. Fig. 2). Therefore, the cosmic-ray fluxes typically lie many orders of magnitude below 
the background, being completely unobservable. 



5 Conclusions 

Cosmological scenarios where the observed matter-antimatter asymmetry is generated 
by the supersymmetric thermal leptogenesis mechanism generically fail to reproduce the 
observed abundances of primordial elements. In the minimal scenario, the LSP must 
be a gravitino heavier than ~ 5 GeV. Therefore, if R-parity is conserved, the NLSP is 
very long lived, jeopardizing the successful predictions of the standard BBN scenario. To 
solve this conflict we have postulated the existence of a light hidden sector fermion which 
couples very weakly to the NLSP. If the coupling is large enough, the NLSP will decay 
dominantly into hidden sector fermions before the epoch of primordial nucleosynthesis, 
avoiding all the nucleosynthesis constraints altogether. We have analyzed the constraints 
on this coupling from dark matter overproduction and from structure formation and we 
have found a wide window of parameters where the cosmological history of the Universe 
can be consistent with baryogenesis through thermal leptogenesis and with BBN. We have 
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presented two concrete models to illustrate the viability of the above mentioned scenario. 
Furthermore, we have discussed some experimental signatures at particle colliders which 
can provide evidence for our mechanism and distinguish it from other scenarios. Examples 
of these are the number of Higgs events and the lepton flavour violation in late stau 
decays. Finally, we demonstrated that the scenario is easily compatible with bounds 
from cosmic-ray observations. 
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Note added 

The authors of Ref . [H] , which appeared simultaneously to this one, independently con- 
sidered a similar mechanism to accommodate BBN bounds and thermal leptogenesis for 
a specific situation where a stau LOSP can decay to a goldstino. 
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